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k> \ The complex-time method for quantum tunneling is studied. In one- 

dimensional quantum mechanics, we construct a reduction formula for a 
Green function in the number of turning points based on the WKB approx- 
imation. This formula yields a series, which can be interpreted as a sum 
over the complex-time paths. The weights of the paths are determined. 
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The imaginary-time method is quite useful for study of the tunneling phe- 
nomena in quantum field theories. Existence of the solutions of the equation of 
motion, instanton and bounce, allows one to apply the semi-classical approxima- 
tion for path-integral. It is known to lead to valid results, which one can confirm in 
the WKB approximation in quantum mechanical models. Its range of applicability 
is, however, somewhat limited: In the case when there are degenerate vacua, the 
ordinary instanton calculation gives the right level-splitting only for the lowest- 
lying states. In the case of the decay of the metastable states, the decay formula 
is valid for the lowest meta-stable state, the false "vacuum". These points were 
investigated by many authors in recent years, especially in relation with tunneling 
phenomena in high energy scattering, namely, the baryon and lepton number viola- 
tion. As a result, interacting-instanton method and valley methods ° became 
known to be useful. 

Rubakov, et. al. proposed use of a complex-time method as an alternative. 
In such a formalism, the external lines "live" in real-time, while the tunneling 
is carried out in imaginary time, thus the interplay between the initial conditions 
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and the tunneling configuration is expected to be straightforward, unlike the cases 
when the time is purely imaginary. 

Complex-time method is studied extensively in toy models in the one- 
dimensional quantum mechanics. It was argued to overcome the many short- 
comings of the imaginary time method. However, even in the simple toy models, 
the formalism is full of the ambiguities and mysteries. It is the purpose of this letter 
to expose those points and solve those mysteries, thereby establishing the sound 
basis of the complex-time method, on which decent treatment of the quantum field 
theory should be done. 

In its most successful formulation, '' one considers the analytic continuation 
of the time-integral in Fourier-transform of the the fixed-energy Green functions 
(resolvents). We consider the one- dimensional quantum mechanics of a particle 
of unit mass in a potential V(x) assumed to be smooth enough to allow WKB 



approximation. We assume that it has asymptotic regions I and II and also in a 
wide intermediate region III as in Fig.l. The retarded resolvent (the fixed-energy 
Green function) is defined by the following; 

CO 

G R (x u x f , E) = ( Xi | E+ - 6 _ H \xi) = -iJ dTe^W* (xt | e~ iHT \ Xi ) . (1) 

o 

This is useful as energy spectra is obtained from the poles and cuts of the above. 
Complex-time is introduced by deforming the T-integration contour in (1): If any 
saddle point (s) exist in the complex T-plain, one tries to deform the integration 
contour so that it goes through the saddle point (s). 

Easiest to find among the saddle points are the physical ones: Using the path- 
integral expression 

x(T)=x { 

(x{\e- iHT \xi)=i f Dxe iS , (2) 

The resolvent (1) is written as a double-integral over T and x(t). One then re- 
quires saddle-point condition for both T and x(t). The former leads to the energy 
conservation law and the latter the equation of motion. If the initial point X[ is 
in the allowed region, a solution x(t) starts from t — and moves along the real 
axis of t. When it hits a turning point, it can turn around or start moving into the 
forbidden region with t moving along the imaginary axis. Corresponding to the 
choice of the move at turning points, there are infinite of such saddle points. It 
was proposed that all such saddle points should be taken into account with specific 
weights. 

Closer look of this method reveals many unanswered questions. The complex 
T-plane is plagued with singularities and infinite number of saddle points. Among 
saddle points, there are the ones that cannot be obtained from the above proce- 
dure (which we shall call "unphysical saddle points"). Both of the physical and 
unphysical saddle points are distributed on infinite lattice structures. How the 



path is deformed to avoid the singularities and go through only the physical sad- 
dle points is unknown. Instead, one simply assumes that all the physical saddle 
points contribute. Even so, the weight of each saddle point is a mystery. They are 
determined so that results agree with that of that of the WKB approximation. It 
is not known how they come about from the path. 

In order to investigate the situation, we shall construct the Green function from 
the complete orthonormal set of the WKB-eigenfunctions of the hamiltonian. For 
a state with the eigenvalue A in the continuum spectrum, the second-order WKB 
approximation yields the following; 

A x e { C p{x '^ dx ' + B x e-' 1 C PWA for x E I, 



4> x (x) = . x { (3) 

V2^) I ( Cx e- l ^ p{x ' )dx ' + D x eW )d *) for x G II, 





where p(x) = y 2(A — V(x)). Among the coefficients A, B, C, and D, only two are 
independent. Their inter-relation can be summarized in the following relation; 



(4) 



The 2x2 matrix 5(A) is determined by the shape of the potential V(x) in the inter- 
mediate region. The flux conservation law allows the following general parametriza- 
tion; 

(e ia cosh p e % P sinh p \ 
(5) 
e l P sinh p e ta cosh p J 

where a, (3 and p are real functions of A. For a given value of A, there are two 
eigenf unctions. We can show that one orthonormal choice is the following; 

— — , 0, 1, -e^-fltanhp , 
cosnp / 



(ia \ 
e^^)tanhp, 1, 0, —^-) 



(6) 



The Green function is written as an integration over continuous A and sum over 
the discrete spectrum. This results in the following expression in case X[,X{ Gil as 
in Fig.l; 

G(x u xf,E) = -b(zi)p(sf)r 1/2 e- Ai (e Af + iRe~ A{ ) . (7) 

where the argument is A = E and R = iS 2 \/S 22 is essentially the (analytically- 
continued) reflection amplitude for the incoming wave from the right asymptotic 
region II. The exponents are given by 

•X'i.f 

A i)f = J dx\p(x)\. (8) 

Due to the existence of the intermediate region III, the matrix S, which con- 
nects the regions I and II, can be written in terms of the matrices Si that connects 
regions I and III and S 2 for III and II. More specifically, we apply the linear WKB 
connection formula for the latter region and obtain 

fU- Al cosW 2 e- Al sinW 2 \ 
S=S 1 [ 2 . . (9) 

\-e Al smW 2 2e Al cosW 2 ) 

where W 2 and Ai are defined by, 

&2 0,2 

w^JtMx), a^/^!pW|. do) 

a-2 hi 

Using the above, we arrive at the following expressions; 

1 _ 1d iP -2Ai • 

-iR 2 = -i ff* 2Ai = H) + m)e'^ + l-mfe-^ + ..., (12) 

1 + ^R\e ZiXl 2, 



where R\ is the reflection amplitude at the turning point b 



We prove that the reduction formula (11) and (12) is equivalent to summation 
over the classical complex-time paths: Let us first look at the expansions of the 
resolvents in (11). [The convergence of this expansion is guaranteed by 5.] The 
phase W2 is equal to ET + S of the real-time solution between the turning points 
d2 and &2- The factor i/2 is the phase and weight at the turning point, where 
the path approaches from the forbidden region and turns around. When the the 
path hits the turning point from the allowed region and turns around, it picks 
up the phase —%. Thus the series can be understood as the sum over the paths 
that goes through the allowed region (02,62) different times (see Fig. 2). Similarly, 
Ai is — Et + Se and the expansion in (12) can be understood as the summation 
over paths for the forbidden region (61,02) (Fig.3). The reflection coefficient R\ 
itself can be expressed as the sum over the oscillation in (01,61). Combined, they 
produce the sum over the different ways the paths goes back and forth between 
(01, 62). We thus prove that the sum is only over the physical paths. At the same 
time, the weights and phases are completely determined by these expansions. We 
have also looked at the case when X[ and Xf belong to different asymptotic regions. 
We have proven also that it can be obtained by the sum over the physical paths 
with the specified formula for the phase and the weight. 

As a summary, we have shown the Green function is cast in the form of sum only 
over the physical paths. We have derived the general rule for the weights and phases 
for each paths. These coincide the rules assumed by Carlitz and Nicole '' and 
correctly reproduce the WKB energy-level condition and the decay formula. As 
the above procedure is reductive in its nature, our proof is valid for potentials with 
more turning points. Although we have used the linear WKB connection formula 
in this letter, using the quadratic formula, for the purpose of investigating the 
lower energy spectrum, should be straightforward, as our method is solely based 
on the reduction of the matrix S. We should stress that the weight 1/2 for the 
turn around in the allowed region is problematical for the path-integral: How it is 
obtained from a path that goes through all the saddle points is unknown. What 
we have done here is to justify a hybrid method, where we apply the path-integral 



method for the elementary processes and then combine then using the knowledge of 
the wavefunction. In view of the fact that the weights are now fixed by our analysis, 
the path-integral method should be more seriously considered, as this should open 
ways for decent treatment of the tunneling phenomena in field theories. Detail of 
the analysis will be published in near future. 
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FIGURE CAPTIONS 

1) A potential V(X) with asymptotic region I and II and wide intermediate 
region III, where WKB approximation is valid. 

2) Diagrammatic view of the expansion of (11) in terms of i?2- 

3) Diagrammatic view of the expansion of (12) in terms of R\. 
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